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We derive a closed photo-counting formula, including noise counts and a finite quantum efficiency, 
for photon number resolving detectors based on on/off detectors. It applies to detection schemes 
such as array detectors and multiplexing setups. The result renders it possible to compare the 
corresponding measured counting statistics with the true photon number statistics of arbitrary 
quantum states. The photo-counting formula is applied to the discrimination of photon numbers 
of Fock states, squeezed states, and odd coherent states. It is illustrated for coherent states that 
our formula is indispensable for the correct interpretation of quantum effects observed with such 
devices. 
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I. INTRODUCTION 

Quantum optics delivers a manifold of nonclassical fea- 
tures of quantum light. These nonclassical properties can 
be used for various applications in the field of Quantum 
Information Technology [1], Quantum Metrology [2], or 
for highly sensitive measurements, see e.g. [3]. For this 
reason, a reconstruction of at least some properties of 
the quantum state is essential, cf. e.g. [4-6]. One major 
aspect of a quantum state is its photon number distri- 
bution [7, 8]. Sub-Poisson photon statistics renders it 
possible to go beyond the shot noise limit of classical 
states. Even phase sensitive measurements, such as ho- 
modyne detection with a weak local oscillator, require 
the resolution of the photon number distribution [9, 10]. 

However, photon number resolving (PNR) detectors 
are not always accessible. For example, an avalanche 
photo-diode in the Geiger mode is an on/off detector [II]. 
It cannot deliver any information beside the absence or 
presence of photons. Even this tiny bit of information 
can be hidden when taking into account noise and losses, 
see e.g. [12-14]. Surprisingly, a subsequent data analysis 
is possible to obtain a photon statistics with suppressed 
noise and losses [15, 16]. 

There are some proposals for PNR detectors based on 
on/off detectors. In Fig. 1, the general idea is given. The 
incoming light is split into weaker signals. Each of those 
weaker signals has only a small probability to include 
more than one photon. Thus, it can be approximately 
measured by on/ofF detectors. So far the reconstruction 
of the true photon statistics has been done via iterative 
maximum-likelihood methods, see e.g. [17]. The masured 
probabilities can be also directly used to determine, for 
example, photon number correlations [18]. 

One experimentally realized method is the use of detec- 
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FIG. 1: (color online) From the left side quantum light enters 
the measurement scheme. The incident light and A'^ — 1 ports 
with vacuum input are combined by a A'^-port unitary opera- 
tion. The A'^ output beams are detected via on/off detectors. 
The individual clicks are added up to the total number of 
counts. 



tor arrays [8, 19]. Each individual detector in the array is 
an on/off detector. However, simultaneous clicks of some 
detector deliver additional information about the photon 
number distribution of the incident quantum state. An- 
other experimentally realized PNR detector is the multi- 
plex detection scheme, which is based on a multiple 50:50 
splitting of a signal into weaker ones [20-25] . The incom- 
ing signal can be split and detected in different spatial or 
temporal - time-multiplexing - modes or bins. 



In this article we will derive a closed analytical count- 
ing formula for PNR detectors based on on/off detectors. 
In addition, finite quantum efficiencies and noise count 
rates will be included in this new formula. This closed 
counting formula can be used to predict the outcome, 
for example, of a multiplexing or array detector for ar- 
bitrary quasimonochromatic quantum states. We apply 
this method to Fock states, squeezed states, and odd co- 
herent states in different realistic scenarios. We also dis- 
cuss the interpretation of data from PNR detectors for 
determining quantum effects. 

The article is structured as follows. In Sec. II, we de- 
rive the counting formula for on/off detector systems. 
Examples are given in Sec. HI for perfect and imperfect 
measurements. Systematic error effects of on/off detec- 
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tor setups are discussed in Sec. IV. A summary and 
conclusions are given in Sec. V. 



II. PHOTO-COUNTING FORMULA 

Any quantum state can be written in terms of the 
Glauber-Sudarshan P function as [26, 27] 



p = j (faP{a) \a){a\. 



(1) 



Thus, it is sufficient to consider the effects on pure co- 
herent states only. The usual photo-counting formula for 
the true photon statistics reads as 
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where h is the photon number operator, k represents 
the number of photons, and : • ; is the normal order- 
ing prescription, cf. [4, 5]. Taking into account a finite 
quantum efficiency 77 < 1 and noise counts v > Q, the 
photo-counting formula is 
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{rjh + vY 
k\ 
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(3) 



In the following we derive a related closed expression 
for PNR detectors based on on/off' dcitectors. As noted 
above, we restrict our consideration to quasimonochro- 
matic fields, which can be obtained as superpositions of 
initially monochromatic modes. The photo-counting for- 
mula for broad-band radiation is also given in [4, 5]. It 
could be used to further generalize the approach given 
below. 

For this purpose let us consider the incident coherent 
state I a), which is transformed by a unitary operation 

U(Ar), 



aout = U(7V) a;, 



(4) 



with ain/out = (ai,in/out, • • • ,OAr,in/out) being the anni- 
hilation operators for the input/output modes, respec- 
tively. We obtain the input-output relation as 



|a,0, ...,0) I-*- 1^) = \uia,...,UNa), 



(5) 



with {ui)fLi being the first row of the unitary matrix 
U{N), and therefore = 1. In the following we 

consider the detection of this output light field. For a sin- 
gle perfect on/off detector, the projection operator valued 
measure can be given by 

^0 = |0)(0| = :e-": and ^1 = :i - e"":, (6) 

or equivalently 



(7) 



for m = 0(1) for no(one) click, respectively. A set of N 
on/off detectors can be written as 



7r„ = :(e--i)'"^ e"- 



(e"«_i)™~e-"":. 



(8) 



together with the multi-index notion m = (mi, . . . , mjv) 
and rrii G {0, 1}. An output signal of k clicks would mean 
that k detectors deliver a click, |m| = m, = k, 

N 

n. = : E 0e--(e--i)'"^. (9) 

|m|=fe i=l 

Now, we obtain the measured counting probability for 
the coherent state as 



Pkia) =(V|nfc|V) 

= E 

\m\=k 



n(. 

i=l 



(10) 



In addition we may assume that the incident light is 
split into N modes with uniformly distributed intensi- 
ties, |ztip = 1/A^. This assumption yields 



p,ia) = E (el«l^/^ - l) 



|m| 



|m|=fc 
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For the multiplexing measurement the above assumption 
is justified due to the 50 : 50 beam splitters, which divide 
the light beam in each step. This yields that each of the 
2^ = N bins - s is the number of divisions - has the same 
intensity. For the detector array, a Kohler illumination, 
as it is known from classical optics, guaranties that the 
intensity for each of the d x d = N detectors is almost 
identical. In fact, the studied two types of PNR detectors 
differ only in the practically available number of bins. 

Combining Eqs. (1) and (11), we obtain the counting 
statistics for arbitrary quantum states in a closed form 



Pk= j d^aP{a)pk{a 



N\ / ft \ /- ft \ 
=<^^!(iV3^(^"^) (^-^"-) 

where the latter part of the formula expresses the count- 
ing statistics in terms of normally ordered expectation 
values for k counts (fc = 0, . . . , A^). Note that this statis- 
tics represents a quantum version of the binomial distri- 
bution. 

For the inclusion of imperfections, we ristrict our at- 
tention to the quantum efficiency and noise counts, such 
as dark counts or scattered light. For the following we as- 
sume that a given and previously characterized detector 
is sufficiently well described by these disturbances. De- 
pending on technical details of the detector, there may 
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occur other imperfections which could be included if de- 
sired. 

For the presence of losses and noise counts, we derive 
in a similar way 
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(13) 



For this purpose we start with a single imperfect on/ofF 
detector as tt^ = : (e''"+'^ - i)'" 6"''"+'^: instead of 
Eq. (7). The formula (13) is our main finding, represent- 
ing the true photo-counting statistics of multiple on/ofF 
detectors. It delivers the probability for fc counts from 
PNR detectors of the type under study. This expression 
replaces the photo-counting formula (3) for the consid- 
ered PNR detection scheme. Note that, our closed an- 
alytical expression makes numerical procedures dispens- 
able. So far, the latter were the only available meth- 
ods for dealing with such counting statistics, cf. e.g. 
Ref. [12, 13]. 
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FIG. 2: (color online) In each plot, we compare the count- 
ing statistics of Fock states |n), for n = 8 (blue) and 
n = 9 (red) photons. The number of on/off detectors is 
N =10(a), 102(b), 10^(c),10*(d). 



III. EXAMPLES 

In this section, we apply our photo-counting formula. 
First, we derive the Fock basis expansion of Eq. (12). Us- 
ing this new expansion, we consider the discrimination 
of Fock states by N on/off detectors. We compare the 
click-counting statistics with the photo-counting statis- 
tics without noise and losses for a squeezed state. More- 
over, for an odd coherent state we also deal with the 
effects of noise and losses. 



A. Example without imperfections 

First, we consider the discrimination of photon number 
states In this paragraph we aim to study effects of 
the PNR detector itself. Therefore, we may assume for 
the time being a quantum efficiency of 77 = 1 and no noise 
counts, u ^ Q. After some algebra, see Appendix A, it 
can be seen that the projection operator valued measure 
for a PNR detector can be written as 



n 



iV! 

fc!(iV-fc)! 
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(14) 



using the notion in Eq. (11). In Fig. 2, it can be seen 
that the number of on/ofF detectors N must exceed the 
photon number by a few orders of magnitude in order to 
sufficiently discriminate different Fock states. 



Second, let us consider the counting statistics of a 
squeezed- vacuum state |^) [28, 29], 



1 ^/ tenh^ 



|2n). (15) 



This nonclassical quantum state could be the input of a 
multiplexing measurement with iV = 2* bins. In Fig. 3, 
the outcome of a such a measurement for the squeezed- 
vacuum state is shown. It can be seen that the probabil- 
ity to have an odd number of clicks is non-zero. The fact 
that even photon numbers contribute to odd counting 
numbers can be understood as a result of the summation 
of chcks from on/off detectors, cf. Eq. (14). 



B. Example with imperfections 

As an example for the effects of noise and losses, we 
consider the so-called odd coherent state [30] 



l«-)=AA(l«> 



(16) 



-1/2 



with a normalization Af = (2(1 — exp [— 2|q;|^] ) 
The Fock basis expansion of this state only contains odd 
photon numbers, which is opposite to the case of the 
squeezed state. In Fig. 4, we compare the true pho- 
ton statistics based on the Mandel formula (3) with 
the counting statistics from PNR detectors based on 
on/off detectors, cf. Eq. (13). 

It can be seen that the perfect detection with N 
on/off detectors delivers a non-zero probability for mea- 
suring an even click number in the present case, cf. the 
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FIG. 3: (color online) The first figure (a) shows the true pho- 
ton number statistics (green) of a squeezed state with ^ = 1. 
The counting statistics obtained by a multiplexing PNR de- 
tector (blue) for the same state is given for s =3(b), 4(c), 5(d) 
multiplexing steps. Note that the po probability is cut. 



first row of Fig. 4. In this context we have to exclude 
the case of zero clicks, which has still a zero probability. 
This can be easily understood on the basis of Eq. (14), 
which yields Ho — |0)(0| for fc = 0. That is, no clicks can 
only arise from no photons. 



We have also considered the case of imperfections in 
Fig. 4. In the second row we only deal with losses due 
to an imperfect detection efficiency i]. We see that the 
change of both statistics is not much different. This can 
be basically explained by the idea that losses annihilate 
some photons, which in particular affects the statistics for 
small number of counts. Hence it makes not a big differ- 
ence, wether the detector is sensitive to a finite number 
of fc = 0, . . . , clicks, compared with an infinite number 
of fc = 0, 1, . . . photo counts. 

In the third row we deal with the case of noise effects 
only. Here, we find that the change from the perfect to 
the noise case is also qualitatively similar in both cases. 
For sufficiently weak stray light, the noise effects essen- 
tially contribute to the statistiscs at small numbers. Here 
we have to note that a fair comparison of both statistics 
requires that the noise counts of each on/off detector are 
only I'on/off — i^/N. Since we have to add up the signal 
of the PNR detector system over the N detectors, the 
total number of noise counts adds up ly. 

In the fourths row we combine both effects of imperfec- 
tions. Also the combined effects of losses and noise counts 
are qualitatively similar in the two scenarios. In practice, 
the experimenters try to keep dark and stray-light counts 
as small as possible, so they also contribute mainly to the 
events at a small number of counts. Hence the cutoff at 
a finite but sufficiently large number of on/off detectors 
is not very important for the effects of the considered 
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FIG. 4: (color online) The plots (al-a4) show the photon 
statistics (green) of an odd coherent state. In comparison, 
the plots (bl-b4) show the counting statistiscs (red) of PNR 
detectors based on = 25 on/oflf detectors. In the first 
row - figures (al) and (bl) - we consider a perfect detection. 
In the second row, we include only losses -q = 0.8 for both 
detectors (a2) and (b2). In the third row, we consider only 
noise. For having a comparable result, we considered for (a3) 
a noise count u = 0.2 and for (b3) f = 0.2/iV = 0.008. The 
figures (a4) and (b4) include both kinds of the considered 
imperfections. 



imperfections. From these arguments it becomes clear 
that the main differences of both detection scenarios are 
expected to occur due to the cutoff relative to the signal 
field to be detected, rather than due to the considered 
imperfections. 



IV. SYSTEMATIC ERROR EFFECTS 

Let us address the question of systematical errors from 
PNR detectors. For example, the true photon statistics 
of a quantum state can be used to determine nonclassi- 
cal effects. A boundary for classical states can be formu- 
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lated, 
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1 >0, 



(17) 



which is known as the Mandel Q parameter [31]. The 
violation of this classical condition, Q < 0, refers to as 
sub-Poisson photon statistics. However, determining the 
Q value from PNR counting statistics delivers mislead- 
ing results, see Fig. 5. The plotted range is given for 
a maximal value of = 1024, which corresponds to a 
multiplexing detector with s = 10, or a detector array 
with d — 32. The true photon number distribution of 
the considered coherent state is a Poisson statistics, with 
(5 = 0. The plotted negativities are a result of the on/off 
detectors, and they must not be interpreted as nonclas- 
sical effects. Note that such effects from on/off detectors 
also appear in measurements of Bell parameters [32] . 




FIG. 5: (color online) The plot shows the Mandel Q param- 
eter for the counting statistics measured by a PNR detector 
with N on/oS detectors. The computed example is a classical 
coherent state with jaj^ — 20. Note that A*' = 1 denotes a 
measurement with a single on/off detector, and its Q value is 
close to -1. 



The Q parameter for a coherent state with a mean 
photon number |a|^ measured by N on/off detectors is 
Q = exp(-|a|2/A;') - 1. The PNR Q parameter tends to 
the expected value zero for A — >■ cxi, which is equivalent 
to the fact that the binomial distribution converges to 
the Poisson one, see Appendix B. From Fig. 5 it follows 
that a number of A^ = 200 bins is not sufficient to infer 
Poisson photon statistics, Q ~ —0.1, from the measured 
PNR statistics for |q;|^ = 20. Remarkably, in this case 



the number of on/off detectors exceeds the mean photon 
number of the state by one order of magnitude. A care- 
ful interpretation of quantum effects requires that such 
systematic errors have to be taken into account, when 
measuring with PNR detectors. The photo-counting for- 
mula for PNR detectors reveals these systematical errors. 



V. SUMMARY AND CONCLUSIONS 

We derived a closed photo-counting formula for pho- 
ton number resolving detectors based on on/off detectors. 
The cases of perfect and imperfect detection have been 
considered, the latter include finite quantum efficiencies 
and noise counts. It has been shown that the resulting 
photo-counting formula can be interpreted as a quantum 
version of the binomial distribution. For large numbers of 
on/off detectors, this distribution converges to the well- 
known photo-counting formula by Mandel. 

The derived photo-counting formula has been applied 
to photon number states, squeezed states, and odd coher- 
ent states. In particular, we considered the discrimina- 
tion of Fock states in dependence on the number of on/off 
detectors. The true photon statistics of a squeezed- 
vacuum state is compared with that from a photon num- 
ber resolving detector. The photo-counting statistics in 
such a case conceals some properties of the photon num- 
ber distribution of the squeezed-vacuum state. The ef- 
fects of noise and losses for on/off detection setups in 
relation to the Mandel formula have been discussed for 
an odd coherent state. 

For the coherent state, we studied the behavior of the 
Mandel Q parameter and its possible interpretations. We 
showed that the resulting counting statistics may infer 
sub-Poisson light, which is caused by the structure of the 
detector system itself. Such systematic errors can be well 
understood on the basis of our counting formula. 
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Appendix A: Computation of POVM elements in Fock basis 



We aim to derive the Fock expansion for the projection operator valued measure of PNR detectors as given in the 
article, Eq. (14), by starting with 
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Appendix B: Limit for N oo 

Let us consider the Hmit for the photo counting formula of a PNR detector. As it is known from classical probability 
theory, the binominial distribution converges to the Poisson one. Thus, the following proof for the quantum analogue 
is almost identical to its classical version. We obtain from the Taylor expansion of Eq. (13): 
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